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PACS 45.70. -n - Granular systems, classical mechanics of 

Abstract. - In particulate systems with short-range interactions, such as granular matter or simple fluids, 
local structure plays a pivotal role in determining the macroscopic physical properties. Here, we analyse local 
structure metrics derived from the Voronoi diagram of conflgurations of oblate ellipsoids, for various aspect 
ratios a and global volume fractions 0^. We focus on jammed static conflgurations of frictional ellipsoids, 
obtained by tomographic imaging and by discrete element method simulations. In particular, we consider the 
local packing fraction 0/, deflned as the particle’s volume divided by its Voronoi cell volume. We And that the 
probability P(0/) for a Voronoi cell to have a given local packing fraction shows the same scaling behaviour 
as function of 0^ as observed for random sphere packs. Surprisingly, this scaling behaviour is further found 
to be independent of the particle aspect ratio. By contrast, the typical Voronoi cell shape, quantifled by the 
Minkowski tensor anisotropy index p - points towards a signiflcant difference between random packings 
of spheres and those of oblate ellipsoids. While the average cell shape p of all cells with a given value of 
0/ is very similar in dense and loose jammed sphere packings, the structure of dense and loose ellipsoid 
packings differs substantially such that this does not hold true. This non-universality has implications for our 
understanding of jamming of aspherical particles. 


The universality of many features of disordered packings of 
spherical beads, with respect to preparation protocols and sys¬ 
tem parameters, is manifest in various properties, such as the 
universal value of the random close packing limit |[Tj and the 
universal distributions for contact numbers free volumes 
and Voronoi cell shape measures While ellip¬ 
soidal particles and other aspherical particles |T 3 }j 3 T| 

are receiving increasing attention, these questions of universal¬ 
ity, including independence of system parameters and prepara¬ 
tion protocols, have not been comprehensively addressed yet. 
A qualitative difference between ellipsoid and sphere packings 
is revealed by the analysis of the Voronoi diagram of ellipsoid 
packings from various experimental and simulated origins. 

Preparation protocols. - The experimental datasets (rep¬ 
resented by the symbol ■) comprise packings prepared by dif- 


o 

Fig. 1: The Set Voronoi Diagram of ellipsoidal particles has, in gen¬ 
eral, non-convex cells and curved facets and edges, {left) Subset of an 
ellipsoid packing with Set Voronoi cells, {centre) Single particle with 
its cell, {right) Eigensystem of the cell. 


ferent protocols (fluidised beds, different funnels, grids, pour¬ 
ing particles, etc.) and compaction by vertical tapping. The el¬ 
lipsoids have half-axes a \ c : c with a < c and the aspect ratio 
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Data set 

symbol 

jammed 

gravity 

friction (p) 

#particles 

Experiment 

□ 

yes 

yes 

0.38 - 0.75 

^ 5000 

Discrete Element Method (DEM) 

A 

yes 

yes 

no 

9323 

Discrete Element Method (DEM) 

♦ 

yes 

yes 

0.01 - 1 

9323 

Discrete Element Method (DEM) 

T 

yes 

yes 

1000 

9323 

Molecular Dynamics (MD) 

o 

no 

no 

no 

512 

Monte Carlo (MC) 

• 

no 

no 

no 

512 

Expansion technique 

€ 

yes 

no 

no 

1025 


Table 1: Properties of the data sets analysed in the present work. See main text for details. 


is defined SiS a = ajc. These datasets are the same as those used 
in ref. comprising in total 73 datasets of jammed oblate 
ellipsoids of 5 different aspect ratios a and two different parti¬ 
cle types (rough 3D printed particles with a = 0.4,0.6,0.8,1.0 
and smoother pharmaceutical placebo pills with a = 0.59); the 
larger half-axis is c = 3mm (smallest particles) to 4mm (largest 
particles). The standard deviation of the particle volumes is 
2-3 %. The packings were imaged by X-ray tomography; im¬ 
age processing was used to extract particle center coordi¬ 
nates and orientations. The packings consisted of approx. 5000 
particles, of which 600-900 were sufficiently far from cylinder 
walls to be included in the analysis; for the sake of spatial ho¬ 
mogeneity, all packings where radial variations of the packing 
fraction about the mean exceed 0.66% were discarded (as in 
ref. p2|). The particle positions and orientations of all experi¬ 
ments reported here can be downloaded from the Dryad repos¬ 
itory p^ . 

We compare our experimental packings to 120 configura¬ 
tions obtained from discrete element simulations (DEM, (T0|). 
These include datasets of ellipsoids sedimenting into a square 
box in a viscous fluid under the action of gravity, both of fric¬ 
tionless particles (a, RCP) giving highest values of (pg and par¬ 
ticles with very high friction coefficient (t. Sedimented Loose 
Packing, SLP) giving looser packings, and intermediate values 
of friction and viscosity (♦). The RCP estimate of the sed¬ 
imented data is consistent with configurations obtained via a 
particle expansion method (C) 

Reference data for equilibrium configurations of hard-core 
ellipsoids without gravity in the isotropic fluid phase (Fig. 
was obtained by event-driven molecular dynamics simulations 
p5| (o, MD, the same data sets as in ref. pfij ) and by canonical 
Monte Carlo simulations (•, MC, p7|). 

Local structure metrics from the Set Voronoi diagram. - 

For aspherical particles the Set Voronoi diagram p8] - |^ pro¬ 
vides a natural definition of the partition of space into N cells, 
each containing one of the N particles. The Voronoi cell of par¬ 
ticle i is the compact set of points closer to particle i than to any 
other particle. The distance from a point in space to a particle 
is measured as the Euclidean distance to the nearest point on 
the bounding surface of the particle. This is in contrast to the 
conventional Voronoi diagram where the distance is measured 
with respect to the particle centre. Facets of the Set Voronoi 
diagram are in general curved and cells are non-convex, see 
Fig. For monodisperse spheres, the Set Voronoi diagram, 
henceforth simply referred to as the ’Voronoi diagram’, reduces 


to the conventional Voronoi diagram. 

The local packing fraction of particle i is defined as 
((pi)i = yelyi where = 47rac^/3 is the volume of the parti¬ 
cle and Vi the volume of the Voronoi cell Kf containing particle 
i. We characterise the shape of the Voronoi cell Ki by its vol- 
ume moment tensor x<S)xdv where x is the position 

vector relative to the center of mass C/ of K. Similar to the 
tensor of inertia, this tensor captures the distribution of mass; 
the notation derives from the theory of Minkowski tensors 
and integral geometry pTp^ . The three eigenvalues of this 
tensor are The ratio of minimal to maxi¬ 
mal eigenvalue Pi = ^ (0^ 1] is an indicator of the shape 

anisotropy of the Voronoi cell K of particle i. Small values of 
Pi indicate elongated (anisotropic) cells. Note the difference 
to measures of asphericity 0 that quantify deviations from a 
spherical shape; the measure Pi is 1 (and K said to be isotropic) 
for any shape that has statistically identical mass distribution in 
any set of three orthogonal directions; this includes the sphere, 
but also regular polyhedra and the FCC, BCC and HCP Voronoi 
cells (441 . 


Probability distribution of Voronoi cell volumes. - The 

distribution of the Voronoi cell volumes of sphere packs has 
been studied in the context of granular materials ||^[^, super¬ 
cooled liquids |[^, and also with respect to granular entropy and 
the Edwards ensemble p3]^55| |. Aste et al. (^ have shown that 
in random jammed sphere packings below the RCP limit, the 
distribution of Voronoi volumes is universal and independent of 
the preparation protocol. Starr et al 0 have obtained a similar 
result for super-cooled liquids. In both cases, the distributions 
of Voronoi cell volumes collapse when plotted as function of 
(y - <y))/(T where a is the standard deviation of the distribution 
P(y) and (y) its average. Aste et al (^ proposed a derivation for 
a scaling P(fi“^™«)/«y>-y^,„)) where Vmin is the smallest possible 
Voronoi cell of equal-sized spheres and (y) = Yli=i yilN the 
average over all Voronoi cells. 

Figure demonstrates that this universality is not restricted 
to sphere packings, but holds more generally for jammed el¬ 
lipsoid packings: To the resolution of our data, the functional 
form of the distribution depends neither on the global packing 
fraction (pg nor on the particle aspect ratio a. The probability 
for a Voronoi cell in a jammed configuration with global pack¬ 
ing fraction cpg to have local packing fraction (pi is written as 
P{(pi I (pg). When plotted as crP^cpi \ cpg) vs. it is in¬ 

variant for all values of cpg and a, see also Ref. (3^ . This plot 
shows good agreement between the experimental packings and 
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(^,-^g)/G(<Dg) 


Fig. 2: Probability distribution of local packing fractions 0/ of 
jammed static ellipsoid packings (see Table [^. Both the cumulative 
distribution function (main plot) and the probability density (left in¬ 
set), comprising data sets between SLP and RCP with different aspect 
ratios 0.3 < a < 1.0 and friction coefficients 0 < ju < 1000, sup¬ 
port the universality of the probability distribution of Set Voronoi cell 
volumes. The right inset shows the collapse of the standard deviations. 


equilibrium jammed crystal 



global packing fraction (l)g 

Fig. 3: Average anisotropy (J3) of the Set Voronoi cells of the ellip¬ 
soids as function of 0^ for equilibrium ellipsoid configurations, static 
jammed ellipsoid and crystal packings. Dashes on the right-hand ver¬ 
tical axis mark the anisotropies of the particles themselves, i.e. p eval¬ 
uated for a particle rather than its Voronoi cell. 


jammed packings from simulations across the range of accessi¬ 
ble packing fractions (those between SLP and RCP, see Fig.[^ 
and aspect ratios 0.3 < a < 1. By contrast, data from equilib¬ 
rium configurations does not rescale to the same curve. 

The shape and anisotropy of the typical Voronoi cell 
(global averages). - Treating the Voronoi cell volume, or 
equivalently the local packing fraction 0/, as the leading term 
of a shape description of the Voronoi cells, we now proceed 
to higher-order terms. While other scalar quantities, such as 
surface area, integrated curvatures or asphericities may contain 
signatures of such higher-order terms, we here use the tenso- 
rial shape measure similar to the tensor of inertia, and its 
eigenvalue ratio p to quantify the elongation of a cell. 

Figure shows the average Voronoi cell shape anisotropy, 
quantified by {{3) = ^ as function of global pack¬ 

ing fraction 0^. Data is for equilibrium ellipsoid fluids, exper¬ 
iments and simulations of jammed random ellipsoid packings 
and for two dense crystalline configurations (the stretched fee 
obtained by scaling the x coordinate of the fee sphere packing, 
and the densest known structures discussed by Donev et al Q). 

For equilibrium fluids in the limit of vanishing density 
(pg 0, where the typical distance between particles is large 
compared to the particle size, the Voronoi cell shape is indepen¬ 
dent of the particle shape. Consequently, the shape anisotropy 
corresponds to the value p ^ 0.37 of the Poisson point pro¬ 
cess For denser equilibrium fluids |[^[^ the trend of 

the Voronoi shape anisotropy measure {p) can be understood 
by realising that the shape of the Voronoi cells becomes more 
similar to the shape of the particle itself when fg increases (see 
dashes on the right hand vertical axis in Fig.[^ evaluated for an 
ellipsoidal particle itself, rather than its Voronoi cell, the ratio 


h p = (^/cf = a^, see appendix of ref. 1^). When a is small, 
the curve (/5)(0^) hence decreases, while for larger a, (^)(0^) 
increases with fg. 

For the jammed packings, between SLP and RCP, our results 
for spheres (a = 1) are in quantitative agreement with pre¬ 
viously published data j^, with the cells becoming less elon¬ 
gated upon compaction, i.e. yS increases with increasing fg. For 
ellipsoids with smaller value of a, the slope of P((pg) becomes 
smaller and eventually even adopts slightly negative values for 
small a < 0.60. There is an excellent agreement between the 
experimenal packings (□, with different preparation protocols) 
and the numerical data points from DEM simulations (t, ♦, 
a). As previously found l [4^ , sphere configurations exhibit a 
gap in shape anisotropy between the densest equilibrium con¬ 
figuration and the loosest jammed states. For ellipsoids, this 
discontinuity shrinks as the particle’s aspect ratio decreases. 

The shape and anisotropy of the typical Voronoi cell of 
a given size (local analysis). - The packing fraction fg rep¬ 
resents a useful global parameter, which is easily accessible 
in experiments. However, there is no conceivable mechanism 
by which a locally defined quantity, such as the Voronoi cell 
shape or also the contact number, can depend directly (i.e. by 
an immediate causal relation) on 0^; a particle cannot sense the 
global packing fraction. That said, in packings with sufficient 
spatial homogeneity, correlations between averages of the lo¬ 
cal shape metrics and the global packing fraction are evidently 
possible, and commonly observed. Specifically, the study of 
the average contact number Z as a function of fg is a founda¬ 
tion of the jamming paradigm l [3^[57| . 

Here, we use a local density-resolved analysis based on the 
idea that the physical mechanisms underlying granular matter 
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(a) Voronoi diagram with cells colored according to 
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(b) Cell volume distribution P(0/) = P(0/|0^) 



(c) Shape distributions P{P) = P{j3\(j)i^a^ X) 



Fig. 4: The concept of the density-resolved local structure analysis, 
here illustrated w.r.t. (5 for a 2D system. Data in this plot are gener¬ 
ated from 100 configurations of Monte Carlo simulations of a planar 
equilibrium hard-core ellipse fluid, of 1944 particles with a - 0.50 
and 0^ = 0.34. Error bars in (d) represent the standard deviation of the 
distributions in (c). 


occur at the particle scale. This idea was applied to contact 
numbers in Ref. ^7\ and is applied here to Voronoi cell shapes. 
Observed correlations between a local structure metric and the 
local packing fraction 0/ are hence more likely to yield physi¬ 
cal insight than those with the global average (pg. A similar ap¬ 
proach has been used for the analysis of sphere packings @ 0 - 
Figure [^illustrates the concept of the local density-resolved 



local packing fraction (j). 


Fig. 5: Relationship between local packing fraction 0/ and the Set 
Voronoi cell anisotropy p. Each curve represents data for a fixed 
value of (pg, and is averaged over three independent realisations, 
a) Local analysis for spheres (n = 1) and slightly oblate ellipsoids 
(a = 0.8). Gray data points represent experimental configurations, for 
all available values of (pg between RLP and RCP, and the red and blue 
curves show DEM simulation results for SLP and RCP configurations. 
Curves on the bottom (using the right hand abscissa) show P(0/) for 
a = 0.8. b) Local Voronoi cell anisotropy {p){pi) for the DEM esti¬ 
mates for SLP and RCP packings. 


analysis. Particles are grouped by their local packing frac¬ 
tion (0/)/, i.e. into sets <S(0/) composed of all particles i with 
(pi - ^/2 < (0/)/ < 0/ + ^/2 for (pi = A, 2A, 3A,... with a small 
interval A (A = 0.1 in Fig. [^ A = 0.02 in Fig. [^. We define 
the function P{p \ (pi, a, X), which is the probability distribu¬ 
tion of the shape measures p, restricted to the cells in <S(0/), i.e. 
to those with local packing fraction 0/. The unknown param¬ 
eters X capture influences from the packing protocol, friction 
etc. As a result, the X can correlate with (pg even though there 
need not be a causal dependency of the X on (pg. The average 
{P){(pu C(,X) = f p P(p I 0/, a, X) dp over all cells in <S(0/) pro¬ 
vides information on how local structure changes depending on 
local packing fraction 0/. In general, (/?) also depends on the 
aspect ratio a and the unknown parameters X. 

Figure shows the result of this local structure analysis of p 
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of jammed ellipsoid configurations. The key result is the fol¬ 
lowing difference between sphere and ellipsoid packings: in 
sphere packings, the average shape of the Voronoi cells of a 
given local packing fraction 0/ is, as far as it is captured by 
the anisotropy index {3, almost identical in dense and loose 
packings. This is evidenced by the near-collapse of the curves 
a = 1,Z) for packings of different global packing frac¬ 
tion. (j8} is a function of 0/ only, but is largely independent of 
the unknown parameters X, the packing protocol and the parti¬ 
cle friction. 

In ellipsoid packings, illustrated for a = 0.8 in Fig. |^a), 
the curves for different 0^ do not collapse. The average 
{yS)(0/,a,X) depends on both a and X. This indicates that 
packings with low and high 0^ exhibit differences in their lo¬ 
cal structures controlled by a and X. Figure [^b) demonstrates 
the validity of this result for other aspect ratios. It is evident 
that the densest and loosest simulated packings have different 
structures, except around a = 1.0 (spheres) and a = 0.6 (close 
to the densest random ellipsoid packing). 

Discussion and Conclusion. - We have analysed the 
Voronoi diagram of oblate ellipsoid packings, establishing 
which aspects of the Voronoi diagram are universal, i.e. in¬ 
dependent of preparation protocol and particle aspect ratio a 
and further parameters X, and which ones are not. Consider¬ 
ing the geometric nature of this packing problem, these results 
have ramifications for our understanding of jammed systems 
and disordered solids. 

The fact that ellipsoidal particles produce denser random 
packings than spheres is well established, with quantitative 
agreement between different studies of the value of the ’ran¬ 
dom close packing’ limit ^rcp{(^) as function of aspect ratio 
for oblate ellipsoids see the curve labelled RCP in 

Fig. However, while a mean-field theory for (pRCpic^) is de¬ 
veloping GD’ an intuitive geometric understanding for (ppcpict) 
is lacking. In this regard it is noteworthy that at a 0.65, the 
aspect ratio for which (pRCp{(x) is highest | |^10||16| , the Voronoi 
cell shapes are found to be independent of the packing fraction. 
As far as captured by {3 of the volume moment tensor the 
shapes remain approximately constant for all jammed packings, 
both in the global (Fig. and in the density-resolved analysis 
(Fig.igb)). 

The results of Fig. [^emphasise an important distinction be- 
tween random packings of spherical beads and those of aspher- 
ical beads. The structure of spherical bead packs is universal in 
the following sense: on average, the local structure of the typi¬ 
cal particle of a given fixed but arbitrary local packing fraction 
0/ is very similar in differently prepared packings, in partic¬ 
ular with different 0^. This observation, here made w.r.t. the 
Voronoi cell anisotropy of the volume moment tensor, is con¬ 
sistent with similar results for local contact numbers ||^[3^. 

It implies that, at least with respect to averages of the vol¬ 
ume tensor shape measure, the following interpretation of ran¬ 
dom jammed sphere packs is feasible. We consider a pool of 
local structure motifs for each value of the local packing frac¬ 
tion 0/, given by the distributions in Fig.|^c). For spheres (but 
not for ellipsoids), these pools are universal in the sense that. 


for a fixed value of 0/, the same pools can be used to construct 
packings of various global packing fractions 0^. A jammed 
configuration can then be thought of as the composition of ran¬ 
domly drawn elements from the pools; the probability distri¬ 
bution P{(pi) determines the fraction of cells to be drawn from 
each 0/-pool. While this clearly does not represent a construc¬ 
tive approach for the generation of disordered bead packs, it 
illustrates the universal nature of the sphere packing problem, 
manifest in the fact that the same pools of structural elements 
are used for all global packing fractions, just in different pro¬ 
portions. For ellipsoid packings, this universality breaks down 
and motifs in the 0/ pools depend on further parameters X and 
hence correlate with the global packing fraction. We speculate 
that this geometric non-universality is likely to be paralleled by 
a significantly less universal nature of the random close pack¬ 
ing problem in aspherical particles. To get a deeper insight into 
the structure formation of non-spherical particles it is important 
to understand the origin of this non-universality. 

Beyond these specific results for ellipsoidal particles, our 
analysis demonstrates the importance of the correct choice for 
the relevant parameters for the discussion of local structure 
metrics in granular matter. An analysis in terms of the local 
packing fraction 0/, which may in principal directly relate to 
local physical processes, is more meaningful than the conven¬ 
tional analysis in terms of the global packing fraction 0^. 

Acknowledgments. - We thank Weimer Pharma GmbH 
for providing placebo pills, and Rolf Schilling for insightful 
discussions. We acknowledge funding by the German Sci¬ 
ence Foundation (DFG) through the research group ’’Geometry 
and Physics of Spatial Random Systems” under grant SCHR- 
1148/3-2. We thank Francesco Sciortino for the MD data. 

REFERENCES 

[1] Scott G. D., Knight K. R., Bernal J. D. and Mason J., Nature, 
194(1962)956957. 

[2] Aste T, Saadatear M. and Senden T. J., J. Stat. Mech. Theo. 
Exp., 2006 (2006) P07010. 

[3] Starr F. W., Sastry S., Douglas J. F. and Glotzer S. C., Phys. 
Rev. Lett., 89 (2002) 125501. 

[4] Aste T, Matteo T. D., Saadatear M., Senden T. J., Schroter 
M. and Swinney H., EPL, 79 (2007) 24003. 

[5] Schroder-Turk G. E., Mickel W., Schroter M., Delaney 
G. W., Saadatear M., Senden T. J., Mecke K. and Aste T, EPL, 
90 (2010) 34001. 

[6] Donev a., Cisse L, Sachs D., Variano E. A., Stillinger F. H., 
Connelly R., Torquato S. and Chaikin P. M., Science, 303 
(2004) 990. 

[7] Donev A., Stillinger F. H., Chaikin P. M. and Torquato S., 
Phys. Rev. Lett., 92 (2004) 255506. 

[8] Baram R. M. and Lind P. G., Phys. Rev. E, 85 (2012) 041301. 

[9] XiA C., Zhu K., Cao Y., Sun H., Kou B. and Wang Y, Soft Mat¬ 
ter, 10 (2014) 990. 

[10] Delaney G. W., Hilton J. E. and Cleary P. W., Phys. Rev. E, 83 
(2011)051305. 

[11] Mailman M., Schreck C. R, O’ Hern C. S. and Chakraborty B., 
Phys. Rev. Lett., 102 (2009) 255501. 


p-5 





F. M. Schaller et al 


[12] Zeravcic Z., Xu N., Liu A. J., Nagel S. R. and Saarloos W. v., 
Europhys. Lett., 87 (2009) 26001. 

[13] Baule a. and Makse H. A., Soft Matter, 10 (2014) 4423. 

[14] Baule Adrian, Mari Romain, Bo Lin, Portal Louis and Makse 
Hernn a., Nat Commun, 4 (2013) . 

[15] Wegner S., Stannarius R., Boese A., Rose G., Szabo B., Someai 
E. and Borzsonyi T., Soft Matter, 10 (2014) 5157. 

[16] Delaney G. W. and Cleary P. W, EPL, 89 (2010) 34002. 

[17] Delaney G., Weaire D., Hutzler S. and Murphy S., Philos. 
Mag. Lett., 85 (2005) 89. 

[18] Williams S. R. and Philipse A. P, Phys. Rev. E, 67 (2003) 
051301. 

[19] Athanassiadis A. G., Miskin M. Z., Kaplan P, Rodenberg N., 
Lee S. H., Merritt J., Brown E., Amend J., Lipson H. and Jaeger 
H. M., Soft Matter, 10 (2014) 48. 

[20] Baker J. and Kudrolli A., Phys. Rev. E, 82 (2010) 061304. 

[21] Haji-Akbari a., Engel M., Keys A. S., Zheng X., Petschek 
R. G., Paleey-Muhoray P. and Glotzer S. C., Nature, 462 (2009) 
773. 

[22] Jaeger H. M., Soft Matter, 11 (2015) 12. 

[23] Neudecker M., Ulrich S., Herminghaus S. and Schroter M., 
Phys. Rev. Lett., Ill (2013) 028001. 

[24] Smith K. C., Srivastava L, Eisher T. S. and Alam M., Phys. Rev. 
E, 89 (2014) 042203. 

[25] Torquato S. and Jiao Y., Nature, 460 (2009) 876. 

[26] Blouwolee j. and Eraden S., Europhysics Letters, 76 (2006) 
1095. 

[27] Borzsonyi T., Szabo B., Toros G., Wegner S., Torok J., Som¬ 
eai E., Bien T. and Stannarius R., Physical Review Letters, 108 
(2012) 228302. 

[28] Eu Y, Xi Y, Cao Y. and Wang Y, Physical Review E, 85 (2012) 
051311. 

[29] Kyrylyuk a. V., Haar M. A. v. d., Rossi L., Wouterse A. and 
Philipse A. P, Soft Matter, 1 (2011) 1671. 

[30] Wouterse A., Luding S. and Philipse A. P, Granular Matter, 11 
(2009) 169. 

[31] Xiao-Dan Z., Cheng-Jie X., Xiang-Hui X. and Yu-Jie W., Chin. 
Phys. B, 23 (2014) 044501. 

[32] Schaller E. M., Neudecker M., Saadatear M., Delaney G., 
Schroder-Turk G. E. and Schroter M., Phys. Rev. Lett., (2015) 
accepted. 

[33] Schaller E. M., Neudecker M., Saadatear M., Delaney G., 
Mecke K., Schroder-Turk G. E. and Schroter M., AlP Conf. 
Proc., 1542 (2013) 377. 

[34] http://dx.doi.org/10.5061/dryad.rf623. 

[35] Michele C. D., Journal of Computational Physics, 229 (2010) 
3276. 

[36] De Michele C., Schilling R. and Sciortino R, Phys. Rev. Lett., 
98 (2007) 265702. 

[37] Peleiderer P, Milinkovic K. and Schilling T., EPL, 84 (2008) 
16003. 

[38] Schaller E. M., Kapeer S. C., Evans M. E., Hoeemann M. J., 
Aste T., Saadatear M., Mecke K., Delaney G. W. and 
Schroder-Turk G. E., Philos. Mag., 93 (2013) 3993. 

[39] Luchnikov V., Medvedev N., Oger L. and Troadec J., Phys. Rev. 
E, 59 (1999) 7205. 

[40] Luchnikov V., Gavrilova M., Medvedev N. and Voloshin V., Eu- 
ture Generation Computer Systems, 18 (2002) 673 . 

[41] Schroder-Turk G. E., Mickel W, Kapeer S. C., Klatt M. A., 
Schaller E. M., Hoeemann M. J. R, Kleppmann N., Arm¬ 
strong P, Inayat a.. Hug D., Reichelsdoreer M., Peukert W., 
ScHwiEGER W. and Mecke K., Adv. Mater, 23 (2011) 2535. 


[42] Schroder-Turk G. E., Mickel W., Kapeer S. C., Schaller E. M., 
Breidenbach B., Hug D. and Mecke K., New J. Phys., 15 (2013) 
083028. 

[43] Mickel W, Kapeer S. C., Schrder-Turk G. E. and Mecke K., 
The Journal of Chemical Physics, 138 (2013) . 

[44] Kapeer S. C., Mickel W., Mecke K. and Schroder-Turk G. E., 
Phys. Rev. E, 85 (2012) 030301. 

[45] Edwards S. E. and Oakeshott R. B. S., Physica A, 157 (1989) 
1080. 

[46] Anikeenko a. V., Medvedev N. N. and Aste T., Phys. Rev. E, 11 
(2008) 031101. 

[47] Zhao S.-C. and Schroter M., Soft Matter, 10 (2014) 4208. 

[48] Song C., Wang P. and Makse H. A., Nature, 453 (2008) 629. 

[49] Pica Ciamarra M., Richard P, Schroter M. and Tighe B. P, 
Soft Matter, 8 (2012) 9731. 

[50] Bi D., Henkes S., Daniels K. E. and Chakraborty B., 
arXiv:1404.1854, (2014). 

[51] Puckett J. G. and Daniels K. E., Physical Review Letters, 110 
(2013) 058001. 

[52] Paillusson E. and Erenkel D., Physical Review Letters, 109 
(2012) 208001. 

[53] Wu Y. and Teitel S., Physical Review E, 91 (2015) 022207. 

[54] Asenjo D., Paillusson E. and Erenkel D., Physical Review Let¬ 
ters, 112 (2014)098002. 

[55] Kumar V. S. and Kumaran V, J. Chem. Phys., 123 (2005) 
114501. 

[56] Kapeer S. C., Mickel W, Schaller E. M., Spanner M., Goll C., 
Nogawa T., Ito N., Mecke K. and Schroder-Turk G. E., J. Stat. 
Mech., 2010 (2010) PI 1010:1 (Note that this referenee uses the 
particle center as the origin for the analysis of , rather than 
the center of mass of the Voronoi cell). 

[57] O’Hern C. S., Silbert L. E., Liu A. J. and Nagel S. R., Phys. 
Rev. E, 68 (2003) 011306. 


p-6 



